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TRANSIENT  SIGNAL  DISTORTION  IN  A  MULTIPATH  ENVIRONMENT 


R.L.  Field 
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The  classi f icat ion  of  transient  signals 
is  facilitated  if  propagat iona 1  distor¬ 
tion  is  removed.  This  distortion  can  be 
removed  if  the  ocean  impulse  response 
from  source  to  receiver  is  known.  The 
ocean  impulse  response  is,  however,  very 
sensitive  to  source  location  and  the 
ocean  environment.  It  is  shown  that  for 
a  bottom-mounted  receiver  and  almost 
normal  incidence,  distortion  of  a  244-Hz 
transient  signal  envelope  can  be 
accounted  for  if  the  time  delays  and 
relative  amplitudes  of  the  ocean  impulse 
response  are  accurately  known.  These 
time  delays  and  amplitudes  are  extracted 
from  transient  data  for  the  known  source 
case.  This  work  provides  the  motivation 
for  the  development  of  h  igh-resolut ion, 
time-delay  estimation  algorithms  to  be 
presented  in  a  companion  paper .  These 
algorithms  will  be  tested  on  the  data 
presented  in  this  paper  for  the  unknown 
source  case. 


convolution.  (Since  the  transient  signal 
is  assumed  to  have  a  high  signal- 
to-noise  .ratio,  no  noise  is  considered 
in  the  equations.)  If  the  source  has 
been  localized  and  ttie  environment  is 
accurately  known,  then  an  approximation 
to  g(t>  can  be  calculated  with  an 
appropriate  acoustic  propagation  model.1 
This  model  response  will  be  denoted  by 
9m<t>. 

The  desired  s(t)  can  be  obtained  by 
deconvolving  r(t)  with  the  inverse  of 
the  impulse  response  computed  by  the 
model,  g^'ft).  Convolving  both  sides  of 
equation  1  with  gm1(t) 


s(t)  -  r(t)  »  (t) 


s(t)  *  [ g ( t )  *  9m  <t>l 


s(t)  *  S(t)  , 


INTRODUCTION 

Transient  source  extraction  can  be  formu¬ 
lated  as  a  time  and  spatially  invariant 
deconvolution  problem  if  the  source  has 
been  localized  in  a  well-known  environ¬ 
ment.  If  these  conditions  are  met,  then 
the  received  signal  r(t),  can  be  repre¬ 
sented  by  the  following  convolution 
integral : 


s(u)  g(t  -  u)  du 


-  s(t)  *  g(t)  , 


where  s(t)  is  the  desired  signal  and 
g(t)  is  the  ocean  impulse  response  or 
Green's  function.  The  asterisk  denotes 


where  s(t)  is  an  estimate  of  the  desired 
source  function  and  5 (t)  is  an  approxima¬ 
tion  to  the  delta  function.  If  6  ( t )  is 
to  be  a  reasonable  estimate  to  S(t)  and, 
therefore,  s(t)  is  to  be  a  good  estimate 
of  s(t),  then  the  environmental  param¬ 
eters  used  in  the  propagation  model  and 
the  source  location  must  be  accurately 
known . 

In  general,  neither  the  environment  nor 
the  source  location  w^i  1 1  be  known 
accurately  enough  for  S(t)  to  be  a 
reasonable  approximation  to  S(t).  A 
major  reason  why  equation  2  fails  is  the 
sensitivity  of  the  received  signal  on 
the  time  delays  of  g(t)  The  time  delays 
considered  here  are  due  to  multipaths. 

This  paper  provides  the  motivation  for 
development  of  the  high-resolution, 
time-delay  estimation  algorithms  discus¬ 
sed  in  a  companion  paper.2  These 
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high-resolution  algorithms  use  the 
approximate  source  location  and  environ¬ 
mental  information  contained  in  the 
model  ocean  impulse  response  to  extract 
an  unknown  source  function  from  the 
received  signal.  In  this  work,  the 
transient  data  on  which  these  algorithms 
are  being  tested  are  discussed  and 
(almost)  optimum  time  delays  and 
amplitudes  which  these  algorithms  should 
find  are  determined . 


mounted  on  the  source.  The  source 
transmitted  the  gated  sinusoid  as  the 
ship  drifted  over  the  bottom  hydrophone 
shown  in  figure  1.  The  range  is  esti¬ 
mated  to  be  100  m. 

The  ocean  bottom  is  characterized  by  a 
thin  sediment  layer  over  a  highly  reflec¬ 
ting  basalt  (see  figure  1)  .  The  sediment 
varies  from  0  to  20  m.  For  this  problem, 
a  10-m  sediment  thickness  was  chosen. 


TRANSIENT  DATA/MODEL 
IMEULSfi- RESEQMSE 
»  * 

Transient  data  were  gathered  in  the 
Atlantic  Ocean ,on  a  bottom-mounted  hydro¬ 
phone  in  780  m  of  water.  The  experimen¬ 
tal  geometry  is  shown  in  figure  1.  The 
acoustic  source  was  at  a  depth  of  40  m 
and  transmitted  the  244-Hz  gated 
sinusoid  shown  in  figure  2.  The  source 
signature  was  recorded  from  a  hydrophone 


ttmt  (ski 

Figure  2.  Source  function:  244-Hz  gated  sinusoid. 


The  received  signals  from  transmissions 
3  minutes  apart  are  shown  in  figures  3a 
and-  3b.  The  signals  received  in  this 
environment  were  modeled  with  the  fast 
field  program,  SAFARI.3  A  broadband 
Gaus*sian  pulse'  is  used  to  model  the 
impulse  response,  g^ft).  The  model  pre¬ 
dicts  the  four  paths  shown  in  figure  1. 
Path  1  is  the  direct  path,  path  2  is  the 
reflection  off  the  basalt,  path  3  is 
surface  reflected,  and  path  4  is 
reflected  off  the  basalt  from  another 
surface  reflection.  The  model  response 
corresponding  to  figure  3a  is  shown  in 
figure  4  with  the  four  paths  labeled. 

It  was  found  that  repeated  calculations 
of  gm(t)  were  not  sufficiently  accurate 


Figure  3.  (a)  Received  data,  and  (b)  received  data  (3  minutes 
later). 
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Tima  (sac) 


Figure  4  Broadband  Gaussian-pOlse  propagated  with  SAFARI, 
source  depth  =  40  m.  range  =  0  1  km. 

0 

to  reconstruct  either  of  the  received 
signals  shown  in  figure  3.  Small  changes 
(two  or  three  sample  points,  2000-Hz 
sample  rate)  in  any  of  ttie  three  time 
delays  (figure  4)  produced  large  vari¬ 
ations  in  the  shape  of  the  received 
signal  being  modeled. 


TIME  DELAYS  AND  AMPLITUDES 


In  order  to  determine  if  the  model 
predicts  the  correct  form  of  g(t),  an 
inverse  filter  of  the  source  function 
(figure  2)  is  computed.  Applying  this 


inverse,  s_1(t>,  to  the  received  signals 


of  figures  3a  and  3b 


-1, 


g(t)  -  r ( t )  *  s  *(t) 


-1, 


-g(t)*(s(t)«s  (t)] 


g(t)  *  Sg(t)  , 


(3) 


where  g(t)  is  the  impulse  response 
estimated  from  the  received  signal  and 
5«.(t)  is  the  band  limited  approximation 
to  the  delta  function  for  the  known 
source  s(t).  The  estimated  impulse 
lesponses  of  the  received  signals  shown 
in  figures  3a  and  3b  are  displayed  in 
figures  5a  and  5b,  respectively.  The 
responses  have  been  bandpass  filtered 
from  150  to  400  Hz.  Since  s(t)  is  accu¬ 
rately  known,  g(t)  should  be  a 
reasonable  estimate  of  g<t). 

Assuming  the  model  response,  gm(t) 
(figure  4),  is  correct  in  general  form, 
the  four  events  labeled  in  figures  5a 
and  5b  should  be  the  desired  time  delays 


-0  5 


Figure  5  (a)  Estimated  impulse  response  function;  g(t).  and 
(b)  estimated  impulse  response  (3  minutes  later) 


and  amplitudes.  These  are  extracted  from 
the  two  estimated  responses  and 
displayed  in  figures  6a  and  6b. 


To  determine  if  these  extracted  paths 
reproduce  the  received  signals,  the 
waveforms  of  figures  6a  and  6b  are 
convolved  with  the  source  function, 
s(t),  shown  in  figure  2.  The  envelopes 
of  these  reconstructed  signals  are 
compared  with  the  envelopes  of  the 
received  signals.  Figures  7a  and  7b  show 
the  reconstructed  and  received  signals 
overlaid. 


CONCLUSIONS 

From  figure  7  it  is  seen  that  the 
extracted  time  delays  and  amplitudes  of 
figures  6a  and  6b  give  good  estimates  of 
the  received  signals.  This  also  demon¬ 
strates  that  the  model  predicts  the 
correct  form  of  the  true  impulse 
response,  but  due  to  errors  in  locali¬ 
zation  and  environment,  the  model  could 
not  be  used  directly  to  model  the 
received  signal. 
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Figure  6  (a)  Extracted  impulse  response,  and  (b)  extracted 
impulse  response  (3  minutes  later). 
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